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We study the critical behavior of the short-range p-state Potts spin glass in three and four dimen- 
sions using Monte Carlo simulations. In three dimensions, for p = 3, a finite-size scaling analysis of 
the correlation length shows clear evidence of a transition to a spin-glass phase at T c ~ 0.273 for a 
Gaussian distribution of interactions and T c ~ 0.377 for a ± J (bimodal) distribution. These results 
indicate that the lower critical dimension of the three-state Potts glass is below three. By contrast, 
the correlation length of the ten-state (p — 10) Potts glass in three dimensions remains small even 
at very low temperatures and thus shows no sign of a transition. In four dimensions we find that 
the p = 3 Potts glass with Gaussian interactions has a spin-glass transition at T c ~ 0.536. 
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I. INTRODUCTION 

There are two main motivation for studying the Potts 
glass: First, the model for which the number of Potts 
states p is equal to 3 is a model for orientational glassesA 
(also known as quadrupolar glasses), which are randomly 
diluted molecular crystals where quadrupolar moments 
freeze in random orientations upon lowering the tempera- 
ture due to randomness and competing interactions. This 
is similar to spin glasses where spins are frozen in random 
directions in the spin-glass phase. However, unlike con- 
ventional spin-glass systems^ with Ising or vector spins, 
the Potts model does not have spin inversion symmetry, 
except for the special case with p = 2 in which case the 
Potts model reduces to the Ising model. In mean-field 
theory (i.e., for the infinite-range version of the model) 
it is found— that the three-state Potts glass has a finite 
transition temperature to a spin-glass-like phase. 

Second, the Potts glass could help us understand the 
structural glass transition of supercooled fluids. In the 
mean-field (i.e., infinite-range) limit, the Potts glass with 
p > 4 has a behavior quite different from that of the 
Ising spin glass since two^ transitions occur in the Potts 
case. As temperature is lowered, the first transition is a 
dynamical transition at Td below which the autocorre- 
lation functions of the Potts spins do not decay to zero 
in the long-time limit. The second transition is a static 
transition at T c < T<\ where a static order parameter, 
described by "one-step replica symmetry breaking," ap- 
pears discontinuously. At the mean-field level, it turns 
out that the the equations which describe the dynamics 
of the p-spin Potts glass for T > T<j are almost identi- 
cal to the the mode-coupling equations; see Ref. and 
references therein, which describe structural glasses as 
they are cooled towards the glass transition. Hence Potts 
glasses with p > 4 and structural glasses seem to be con- 
nected, at least in mean-field theory. 

However, it is important to understand to what ex- 
tent these predictions for phase transitions in the infinite- 



range Potts glass are also valid for more realistic short- 
range models. Here we consider this question by per- 
forming Monte Carlo simulations on the three-state and 
ten-state Potts glass models. 

Earlier work on short-range three-state Potts glasses 
in three dimensions have used zero-temperature 
domain-wall renormalization group methods, 7-8 Monte 
Carlo simulationS)2ii2iii and high temperature series 
expansions —2iiii4 The results of the domain-wall and 
Monte Carlo studies have generally indicated that the 
lower critical dimension di is equal to, or close to, 3. 
The series expansion results tend to be better behaved 
in higher dimensions and less reliable in d = 3, but 
Schreider and Reger— have also argued in favor of di = 3 
from an analysis of their series. Overall, though, it 
seems to us that conclusive evidence for the value of the 
lower critical dimension is still lacking. In this paper we 
therefore investigate the lower critical dimension of the 
three-state Potts glass by Monte Carlo simulations of 
this model in three space dimensions using the method 
that has been most successful in elucidating the phase 
transition in Isingi^ as well as vector— spin glasses: a 
finite-size scaling analysis of the finite-size correlation 
length. 

In four space dimensions, Scheucher and Reger— per- 
formed Monte Carlo simulations, finding that the critical 
temperature T c is finite and a value for the order param- 
eter exponent (3 close to zero, as well as a correlation 
length exponent v consistent with 1/2. These are the 
expected^ values at the discontinuous transition in the 
mean-field Potts model for p > 4. It is surprising that 
these should be found in four dimensions, since, as for 
Ising and vector spin glasses^ the upper critical dimen- 
sion is 6. To try to clarify this situation we therefore also 
apply the finite-size scaling analysis of the correlation 
length to the three-state Potts glass in d = 4. 

Earlier work on the ten-state Potts glass in three di- 
mensions by Brangian et alw^^ found that, in contrast 
to the mean-field version of the same model, both the 
static and dynamic transitions are wiped out. Relaxation 
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times follow simple Arrhenius behavior, and judging by 
the lack of size dependence in the results, the correla- 
tion length is presumed to be small, although it was not 
calculated explicitly. In this paper we calculate the cor- 
relation length and find that it is always less than one 
lattice spacing, consistent with the results of Brangian et 
al. 

In Sec. [n] we discuss the model and the numerical 
method that we have used. We show the results of the 
simulations for the three-state model in Sec. 11111 and for 
the ten-state model in Sec. IIVI Section |V] contains a 
summary. 



II. MODEL AND METHOD 

In the Potts glass, the Potts spin on each site i can be 
in one of p different states, rij = 1, 2, . . . ,p. If neighboring 
spins at i and j are in the same state, the energy is —Jij, 
whereas if they are in different states, the energy is zero. 
Thus the Hamiltonian of the model is given by 



(1) 



in which the summation is over nearest-neighbor pairs. 
The sites lie on a hypercubic lattice with N — L d sites, 
and we consider dimensions d = 3 and 4. The interac- 
tions Jij are independent random variables with either a 
Gaussian distribution with mean Jo and standard devia- 
tion unity, 
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or a bimodal (± J) distribution with zero mean, 



V(Jij) = -W ij -l)+S(J ij +l)}. (3) 

It is convenient to rewrite the Potts glass Hamiltonian 
using the "simplex" representation where the p states are 
mapped to the p corners of a hypertetrahedron in (p— 1)- 
dimensional space. The state at each site is represented 
by a (p — l)-dimensional unit vector Sj, which takes one 
of p values, S A satisfying 



S A • S A ' = 
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(A = 1,2,- •• ,p). In the simplex representation, the 
Potts Hamiltonian is similar to a (p — l)-component vec- 
tor spin-glass model, 



U = - ^2 J ij S * ' S J> 
ignoring an additive constant, where 

J ij ~ p %3 ' 



(5) 
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However, an important difference compared with a vector 
spin glass is that, except for p = 2, if S A is one of the 
(discrete) allowed vectors, then its inverse — S A is not 
allowed. 

In mean-field theory, assuming the transition is contin- 
uous, the spin-glass transition temperature is given by£ 
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where z — 2d is the number of nearest neighbors and 
[•••]av denotes an average over disorder. Actually, in 
mean-field theory there is a discontinuity in the order 
parameter for p > 4 and the transition then occurs at a 
higher temperature^ For the case of p = 10 this is found 
to be afeSi T c = 0.28, rather than Eq. J7J, which predicts 
T c = 0.24. Hence, the mean-field transition temperatures 
for the models discussed in this paper are 
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three-state, d = 3, 
three-state, d = 4, 
ten-state, d — 3. 
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The spin-glass order parameter for a Potts spin glass 
with wave vector k, q piy (k), is defined to be 
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where /i and v are components of the spin in the simplex 
representation and "(1)" and "(2)" denote two identi- 
cal copies (replicas) of the system with the same disor- 
der. The wave- vector-dependent spin-glass susceptibility 
Xsc(k) is then given by 



XSG (k) = iv£[<i<r(k)r 



(10) 
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where (• • • ) denotes a thermal average. 

The correlation length, of a system of size L is 
defined 15 by 



2sin(fc min /2) 



Xsg(O) 
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where k m ; n — (2ir/L,0, 0) is the smallest nonzero wave 
vector. Since the ratio of £,l/L is dimensionless, it satis- 
fies the finite-size scaling form 
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where v is the correlation length exponent, T c is the tran- 
sition temperature, and X a scaling function. At the 
transition temperature, £l/ L is independent of L, and so 
this quantity is particularly convenient for locating the 
transition. Once T c has been obtained from the intersec- 
tion point, one determines the exponent v by requiring 
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that the data collapse onto a single curve when plotted 
as a function of L X I V (T - T c ). 

We are also interested in obtaining the critical expo- 
nent 77 which characterizes the power-law decay of the 
correlation function at criticality. The two exponents v 
and 77 completely characterize the critical behavior, since 
the other exponents can be obtained from them using 
scaling relations. 22 From scaling theory, we expect the 
spin-glass susceptibility xsg to vary as 

X SG = L^c(l^[T-T c }). (13) 

To perform the scaling fits systematically, we use the 
method recently proposed in Ref. |U By assuming that 
the scaling function can be approximated in the vicinity 
of the critical temperature by a third-order polynomial 
(with coefficients Cj, z=l— 4) , a nonlinear fit can be per- 
formed on the data for £j, with six parameters (e;, T c , 
and v) using Eq. (|12fl in order to determine the opti- 
mal finite-size scaling. For the xsg data, we have an 
additional parameter 7/ from Eq. 113(1 . so the fit involves 
7 parameters. It is necessary to choose a temperature 
range within the scaling region to do the fits. We limit 
our temperature range to where results from scaling the 
data in (T - T c ) and (/3 - /3 C ), where (3 = 1/T is the in- 
verse temperature, agree within error bars and quote the 
results obtained from the fits with (T — T c ). To estimate 
the errors we used a bootstrap procedure, described in 
Ref. |23[ If there are N samp samples for a given size, we 
generate iVboot = 100 random sets of the samples (with 
replacement) with N samp samples in each set, do the anal- 
ysis for each set, and take the standard deviation of the 
results from the sets to be the error. We emphasize that 
this gives the statistical errors only. In addition, there 
can be systematic errors, from corrections to finite-size 
scaling which are difficult to estimate from simulations 
on only a modest range of sizes, as is possible here. 

Spin-glass simulations are hampered by long equili- 
bration times, so we use parallel tempering, 24-25 which 
has proven useful in speeding up equilibration. In this 
approach, copies of the system with the same interac- 
tions are simulated at several temperatures. To check for 
equilibration, we use the following relation which holds if 
the couplings are drawn from a Gaussian distribution 
with mean Jo an d unit standard deviation: 

= [< W J av - [( (S ninj f )] ay - J T [(S ninj )] av , 

(14) 

where U = (1/N)(7i) is the energy per spin for a given 
disorder realization. Equation (|14|) is obtained 2 ^ by in- 
tegrating by parts the expression for the energy with re- 
spect to the interactions Jij. Note that the square can 
be omitted in the second term on the right of Eq. (| 1 4|) 
(since 5 ninj only takes values 1 and 0). The first term 
on the right is calculated from two copies at the same 
temperature— i.e., [(^J^'^-l^j^. Hence, if Nt 
is the number of temperatures, the total number of copies 
simulated (with the same interactions) is 2Nt- As the 
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FIG. 1: (Color online) An equilibration plot showing the left- 
hand side (upper curve) and right-hand side (lower curve) of 
Eq. JT]}. The dat a are for the three-state Potts glass in d = 3 
for L = 12 at T — 0.248 as a function of equilibration time 
measured in Monte Carlo sweeps, iV sw . An equal number 
of sweeps are also performed for measurement, so the total 
number of sweeps is 2iV sw , respectively. The two sets of data 
merge at long times, as expected, and then remain indepen- 
dent of time. The inset shows the correlation length £l for 
the corresponding number of Monte Carlo sweeps. 

simulation proceeds, the left-hand side (LHS) of Eq. (f 1 4|) 
approaches the equilibrium value from above, while the 
right-hand side (RHS) approaches the (same) equilibrium 
value from below since the spins in the two copies are 
initially uncorrelated. When both quantities are simul- 
taneously plotted as a function of Monte Carlo sweeps, 
the two curves merge when equilibration is reached and 
remain together subsequently^ Hence requiring that 
Eq. ijTil) is satisfied is a useful test for equilibration. 



An example of the equilibration test is shown in Fig. ^ 
for the three-state Potts glass in d = 3 for L = 12 at 
T = 0.248. In the main panel, the upper curve is the 
energy term on the left-hand side of Eq. l|14f) while the 
lower curve corresponds to the right-hand side of Eq. i|14|) 
as a function of equilibration time measured in Monte 
Carlo sweeps (MCS). Note that one sweep consists of 
one Metropolis sweep of single spin-flip moves on each 
copy, followed by a sweep of global moves in which spin 
configurations at neighboring temperatures are swapped. 
Figure n shows the number of sweeps used for measure- 
ments, -ZV SW ; in addition, an equal number of sweeps have 
been used for equilibration, so the total number of sweeps 
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is 2N SW . We see that the two curves merge in this case 
at around 10 7 MCS and stay within error bars after that, 
indicating that equilibration has been reached. The inset 
of Fig. shows the correlation length ^ as a function of 
Monte Carlo sweeps. The data increase until ~ 10 7 MCS 
where they saturate. This indicates that the equilibra- 
tion time for the correlation length is in agreement with 
the equilibration time determined from the requirement 
that the two sides of Eq. $Tty agree. 

Unfortunately, the equilibration test using Eq. 1141) 
only works for a Gaussian bond distribution. For the 
bimodal disorder distribution, we study how the results 
vary when the simulation time is successively increased 
by factors of 2 (logarithmic binning of the data) . We re- 
quire that the last three measurements for all observables 
agree within error bars. 

For the Potts glass with p > 2, an additional issue 
arises, which is not present for the Ising (p = 2) case: At 
low temperatures ferromagnetic correlations develop^ 
due to lack of spin-inversion symmetry, even for a sym- 
metric disorder distribution (Jo = 0). Since we want 
to study the glassy phase, rather than the ferromagnetic 
phase, we can mitigate this undesirable feature by choos- 
ing a negative value for Jo. A choice of Jo = —1 seems 
sufficient for the ten-state Potts glass^S while for the 
three-state Potts glass, the ferromagnetic correlations are 
small, as we shall show, so we set Jo = 0. 
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FIG. 2: (Color online) Finite-size correlation length £l/L vs 
T for the three-state Potts glass with Gaussian couplings in 
three dimensions. The curves intersect at T c ~ 0.273. The in- 
set shows a scaling plot according to Eq. il I ill using T c = 0.273 
and v — 1.30. The solid line is the third-order polynomial 
used in the fit. 



III. THREE-STATE POTTS GLASS 
A. Gaussian disorder (d = 3) 



We first discuss results for the three-state Potts glass 
with Gaussian interactions. The simulation parameters 
are shown in Table. |U As shown in Fig. [2 the data for 
£l/L intersect at T ~ 0.27. The common intersection 
provides good evidence for a spin-glass transition at this 
temperature. It would be desirable to obtain data in 
this temperature range for larger sizes, but unfortunately 
equilibration times are prohibitively long. Interestingly, 
the same transition temperature T c ~ 0.27 was proposed 
earlier by Banavar and Cieplak 7 ^ who studied the free- 
energy sensitivity to changes in boundary conditions at 
finite temperature using very small system sizes L = 2 
and 3. 

To determine T c and v, we analyze the £l data for 
L > 6 using the method discussed in Sec. |n] for the tem- 
perature range 0.20 < T < 0.32 for all L and obtain 

T c = 0.273(5), u = 1.30(7). (15) 

A finite-size scaling analysis of £l /L according to Eq. i|12|) 
is shown in the inset of Fig. [3 for L > 6. The solid line 
represents the third-order polynomial that we have used 



TABLE I: Parameters of the simulation for the three-state 
Potts glass with Gaussian interactions and Jo = in three 
dimensions. A samp is the number of samples. A sw is the 
total number of Monte Carlo sweeps used for measurement for 
each of the 2Nt replicas for a single sample. An equal number 
of sweeps is used for equilibration; hence, the total number 
of MCS performed is 2A SW . T . is the lowest temperature 
simulated, and At is the number of temperatures used in the 
parallel tempering method for each system size L. 



L 


± mm 


N T 


Asamp 


A sw 


4 


0.2 


18 


9356 


80000 


6 


0.2 


18 


3980 


80000 


8 


0.2 


20 


4052 


655360 


12 


0.248 


28 


352 


16777216 



to approximate the scaling function. The data collapse 
well. 



Next we discuss the results for XSG- We find that al- 
lowing T c to vary, in addition to v and rj, the value of T c 
so obtained is rather different from that obtained using 
data for £l/L. Since we believe the value of T c obtained 
from £,l/L — namely, T c = 0.273 — to be our most accu- 
rate estimate, we fix T c to this value and fit rj and v (and 
the polynomial parameters Cj) to Eq. I|13|) . This gives 
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FIG. 3: (Color online) Scaling plot of xsg/L 2 ~ v vs T for the 
three-state Potts glass with Gaussian interactions in three 
dimensions according to Eq. 1131 with T c fixed to 0.273 (de- 
termined from the data for £l/L shown in Fig. [21 - The fit to 
Eq. JT3J| then gives r\ = -0.03 and v = 0.94. 

v = 0.94(3) and 

t) = -0.03(2). (16) 

We show the corresponding plot in Fig. 

If we analyze %SG allowing T c to vary, as well as the ex- 
ponents v and r), we obtain T c = 0.254(23), v = 0.99(8), 
and r\ = —0.18(16). The transition temperature and ex- 
ponents are consistent, within the statistical errors, with 
those found for the same fit with T c fixed to be 0.273. 
We can also analyze the xsg data by fixing both T c and 
v to the values obtained from £,l/L and so obtain just 
77. This gives 77 = —0.06(2) which is consistent with the 
result in Eq. I|16|) . As discussed in Sec. IIII CI the results 
from £l/L are expected to be more accurate than those 
from xsg- 

The three-state Potts model is somewhat similar to the 
XY spin glass in that the spins point along directions in 
a two-dimensional plane. The difference is that in the 
XY model the spins can point in any direction in the 
plane while in the three-state Potts model they can only 
point to one of three equally spaced directions. In units 
where the standard deviation of J[- is unity [see Eq. 
a transition was founds in the XY spin glass at T* Y — 
0.34(2). Since our interactions Jy differ from the J^ by a 
factor of (p—l)/p [see Eq. |JBJ] this corresponds to T* Y — 
0.23(2) in the units used here. We find a slightly larger 
value for T c in the Potts case, which is reasonable since 
fluctuations in the Potts model are presumably reduced 
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FIG. 4: (Color online) Ferromagnetic susceptibility xfm 
against T for the three-state Potts glass with Gaussian cou- 
plings in three dimensions. The upper curves plotted with 
solid symbols correspond to Jo = while the lower curves 
with open symbols correspond to Jo = —1, where Jo is the 
mean of the Gaussian interactions. 



TABLE II: Parameters of the simulation for the three-state 
Potts glass with Gaussian interactions in three dimensions 
used to compare results with Jo = and Jo = — 1 where Jo 
is the mean of the distribution. See Table for details. 
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Jo 
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A mm 


N T 
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18 


1000 


81920 
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18 


1000 


81920 


8 





0.2 
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2400 


655360 
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0.2 


18 
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81920 
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relative to those in the XY model by the constraint on 
the spin directions. 



As mentioned before, the Potts glass with p > 2 devel- 
ops ferromagnetic correlations at low temperatures even 
with Jo = 0. To see the extent of this, we have performed 
additional simulations, with the parameters shown in Ta- 
ble [H] for Jo = —1 as well as Jo = 0. We have calculated 
the ferromagnetic susceptibility, defined by 

XFM = iV[(|m| 2 )] av , (17) 

where m = iV _1 J^i^i, and show the data in Fig. 0] 
For Jq = 0, the ferromagnetic susceptibility grows slowly 
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FIG. 5: (Color online) Graph of £l vs T for the three-state 
Potts glass with Gaussian couplings in three dimensions. The 
solid symbols correspond to Jo = while the open symbols 
correspond to Jo = — 1. 



FIG. 6: (Color online) Graph of £l/L vs T for the three- 
state Potts glass with Gaussian couplings in three dimensions 
where the mean of the Gaussian distribution is Jo = —1. The 
vertical line is drawn at the transition temperature T c = 0.273 
of the Jo = case (see Fig. |5J . 



with temperature, seems to be bounded from above, and 
only displays a small system-size dependence. We there- 
fore think that ferromagnetic correlations should not se- 
riously affect the results presented earlier in this subsec- 
tion. For J = — 1, the absolute value of xfm is smaller, 
but shows stronger size and temperature dependence. It 
seems that xfm initially increases as T is reduced, but 
then is strongly suppressed as spin-glass correlations de- 
velop for T — > T c ~ 0.27. The spin-glass correlation 
length for both Jq — and Jo = — 1 is shown in Fig. 
The slope of the Jo = — 1 data varies nonmonotonically 
with temperature, reflecting the nonmonotonic behavior 
of the ferromagnetic correlations. By contrast, the Jo = 
data vary much more smoothly. We further plot £l / L for 
Jo = — 1 in Fig- EI The observed nonmonotonic behavior 
of the slope suggests that corrections to finite-size scaling 
will be large and that the data in Fig. [21 should be more 
reliable in estimating T c . Nonetheless, the Jo = — 1 data 
do have an intersection point, although at a lower tem- 
perature, T ~ 0.20, than the value of 0.273 for the Jo = 
case. We interpret these data to indicate that T c ~ 0.20 
for the Jo = —1 model (of course the value of T c can 
vary with Jo). However, data on a larger range of sizes 
at lower temperature would be needed to confirm this. 
Unfortunately, equilibration times are very long at such 
low temperatures, so these calculations were not feasible. 
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T 

FIG. 7: (Color online) Graph of £l/L vs T for the three- 
state Potts glass with a bimodal coupling distribution in three 
dimensions. The curves intersect at T c ~ 0.377. The inset 
shows a scaling plot according to Eq. 11211 using T c = 0.377 
and v — 1.18. 
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FIG. 8: (Color online) Scaling plot of xsg according to 
Eq. il l .'-il with T c fixed to 0.377 (determined from the data 
for £l/L shown in Fig. EJl for the three-state Potts glass with 
bimodal (±J) interactions in three dimensions. The fit to 
Eq. JT3J| then gives 77 = 0.02 and v = 0.91. 



TABLE III: Parameters of the simulation for the three-state 
Potts glass with bimodal interactions and Jo = in three 
dimensions. For further details see the caption of Tabic Q 



L 


T ■ 

A mm 


N T 


A'samp 


iVsw 


4 


0.2 


18 


10000 


5120 


6 


0.2 


18 


4971 


40960 


8 


0.2 


20 


2046 


1310720 


12 


0.320 


20 


550 


4194304 



T c = 0.377(5) 



v = 1.18(5). (18) 
(|12[) is shown in the 



The scaling plot according to Eq. 
inset of Fig. 

Similar to the data with Gaussian disorder, the fits for 
Xsg give a different value for T c from that obtained from 
data for Since we argue that the latter is more 

accurate, we have fitted xsg with T c fixed to the value 
0.377. This gives v = 0.91(2) and 



r, = 0.02(2) . 



(19) 



These results are shown in Fig. Q3J 

If we allow T c to fluctuate, the fits to the data for 
Xsg give T c = 0.390(21), v = 0.87(7), and 77 = 0.14(17). 
The error bars are larger than in the analysis when T c 
is fixed, but the results are consistent with the previous 
analysis. Note that for the Gaussian case, T c estimated 
from ^l/L is greater than that from XSG> whereas it is 
lower for the ± J distribution. However, these differences 
are within the statistical errors and so are not significant. 
As discussed in Sec. IIII 01 the results from are ex- 

pected to be more accurate than those from xsg- We 
can also analyze the XSG data by fixing both T c and v 
to the values obtained from and so obtain just 77. 

This gives 77 = 0.03(2) which is consistent with the result 
in Eq. 

Comparing our results in this subsection with those 
above for the Gaussian distribution, we see that in both 
cases we obtain a finite T c , showing that di < 3. This 
is in contrast to the T = domain-wall renormalization 
group calculations of Banavar and CieplaliS who find 
di < 3 for the Gaussian distribution and di > 3 for the 
± J distribution. One possible reason for this difference 
is the larger range of sizes that we are able to treat in 
the present study. 

Furthermore, the values for the exponent v agree 
within the statistical errors and those for 77 almost agree; 
see Eqs. (fT5|) . lfTS)l. JTSJl, and OJ. The small additional 
difference in the the values for 77 can be attributed to 
systematic corrections to scaling. Hence our results are 
consistent with universality. 



B. Bimodal disorder (d = 3) 

We expect that the value of the lower critical dimen- 
sion and other universal quantities such as critical expo- 
nents are independent of the distribution of the interac- 
tions. To verify this we present in this subsection data 
for the three-state Potts glass in three dimensions with a 
bimodal (± J) distribution for comparison with the above 
results for the Gaussian distribution. The simulation pa- 
rameters are shown in Table. IIIII 

Our results for £,l/L are shown in Fig. [7\ The data 
intersect at T c ~ 0.38, indicating a spin-glass transition. 
Fitting the data for £ L for 0.34 < T < 0.44 and L > 6, 



C. Gaussian disorder (d = 4) 



For the three-state Potts glass in four dimensions we 
first show data for xfm with Jo = in Fig. [5] to see 
whether significant ferromagnetic correlations develop. 
The behavior is very similar to the case in d = 3 (see 
Fig. 0J since xfm remains small and relatively size in- 
dependent. This indicates that ferromagnetic couplings, 
while present, are relatively weak thus we set Jo = in 
the simulations of the three-state model in d = 4. The 
simulation parameters are presented in Table llVl 
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FIG. 9: (Color online) Graph of xfm vs T for the three- 
state Potts glass with Gaussian couplings with Jo = in four 
dimensions. 




FIG. 10: (Color online) Correlation length £l/L vs T for 
the three-state Potts glass with Gaussian couplings in four 
dimensions. The curves intersect at T c ~ 0.54. The inset 
shows a scaling plot according to Eq. 11211 using T c = 0.536 
and v = 0.81. 



TABLE IV: Parameters of the simulation for the three-state 
Potts glass with Gaussian interactions and Jo = in four 
dimensions. For further details see the caption of Table [I] 



L 


T ■ 

x mm 


N T 


Asamp 


AW 


3 


0.2 


18 


1044 


40960 


4 


0.2 


18 


2691 


40960 


5 


0.350 


17 


1000 


524288 


6 


0.408 


18 


700 


327680 


8 


0.490 


22 


623 


524288 



Next we plot £,l/L in Fig.^J where the data intersect 
at T c « 0.54. Fitting the data for £ L for 0.49 < T < 0.59 
and L > 4, we obtain the best fit as 



T c = 0.536(3), v = 0.81(2) 



(20) 



The scaling plot according to Eq. (|12fl is shown in the 
inset of Fig. [TUI 

Our value for T c is very different from that of Scheucher 
and Regeuii who estimate T c = 0.25(5). However, in the 
vicinity of their estimate for T c they could equilibrate 
only two small sizes L = 3 and 4. Furthermore, they 
used the Binder ratio^ (a ratio of moments of the order 
parameter) which even fails to intersect at the known 
transition temperature in the mean-field Potts model. 
Hence we feel that our approach, using the correlation 
length, is more reliable, especially since we are able to 
use parallel tempering to simulate larger sizes at lower T 
than was possible for Scheucher and Reger in the "pre- 
parallel tempering days." 

We next discuss the data for xsg- First we fix T c to be 
the value 0.536 obtained from Fitting to Eq. H13|) 

then gives v = 0.70(1) and 



n = -0.08(2). 



(21) 



This is shown in Fig. ^] We can also analyze the xsg 
data by fixing both T c and v to the values obtained from 
£l/L, and so obtain just rj. This gives rj = —0.08(2) 
which agrees with the result in Eq. (|21|) . 



However, if we allow T c to vary, as well as v and 77, we 
obtain T c = 0.505(14), v = 0.75(3), and rj = -0.46(16). 
We see that T c is substantially lower than in Eq. I|20|l . 
The resulting value of v is not very different than that 
obtained in Eq. (|20|l but the value of r\ is significantly 
more negative than that obtained in Eq. 1)21(1. These 
differences are due to corrections to finite-size scaling. 
A related problem of different exponents from xsg and 
£,l/L was also found^ for the Ising spin glass (although 
there the main discrepancy involved v rather than rj). 
They argued that the data for £l/L is the most accu- 
rate way to determine T c and v. Since it is dimension- 
less, the data give clean intersections without having to 



9 




FIG. 11: (Color online) Scaling plot of xsg for the three-state 
Potts glass with Gaussian interactions in four dimensions ac- 
cording to Eq. 113H with T c fixed to the value 0.536, which is 
obtained from the data for £l / L shown in Fig. 1101 The fit to 
Eq. JT3J| then gives rj = -0.08 and v = 0.70. 



FIG. 12: (Color online) Ferromagnetic susceptibility xfm vs T 
for the ten-state Potts glass with Gaussian couplings in three 
dimensions. For Jo = —1, represented by the open symbols, 
Xfm remains small and size independent. In contrast, xfm 
for Jo = (solid symbols) becomes very large and shows 
substantial finite-size effects. 



divide by an unknown power of L. Reference |23j con- 
firmed this supposition by reanalyzing their data in an 
alternative way proposed in Ref. 30. The results from 
£l/L did not change significantly, but the value for v 
from xsg changed considerably and became much closer 
to the value found from £l/L. 

Hence, for the present case, we feel that the values 
of v in Eq. l(2Tj|) (obtained from £l/L) and rj in Eq. it2*T|) 
(obtained from xsg with T c fixed at the value from £l / L) 
are the most accurate ones. An estimate for the order 
parameter exponent (3 is then obtained from the scaling 
relation 



which gives 



= (d - 2 + n) 



P = 0.78(3) , 



2 ' 



(22) 



(23) 



This is very different from the value j3 ~ obtained by 
Scheucher and Reger«il Therefore, we find that the tran- 
sition in the three-state Potts glass in four dimensions 
is actually not in the mean-field universality class^ for 
p > 4. 



IV. TEN-STATE POTTS GLASS IN d = 3 
(GAUSSIAN DISORDER) 



o.4 - 




FIG. 13: (Color online) Correlation length £l (not £,l/L) 
against T for the ten-state Potts glass in three dimensions. 
The main panel shows the case for Jo = — 1 where there are 
hardly any finite-size effects. The inset shows the same for 
Jo = showing the size dependence of £l in this case. 
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TABLE V: Parameters of the simulation for the ten-state 
Potts glass with Gaussian interactions with Jo = —1 in three 
dimensions. The second set of parameters for L = 6 and 8 is 
for additional runs to probe the very low-temperature behav- 
ior of the model. 



L 


T ■ 
± mm 


N T 


Asamp 


iVsw 


4 


0.2 


14 


1000 


81920 


6 


0.2 


17 


1000 


131072 


8 


0.2 


19 


998 


262144 


12 


0.2 


26 


343 


20480 


6 


0.02 


28 


469 


262144 


8 


0.08 


18 


600 


131072 



Finally, we turn our attention to the ten-state Potts 
glass in three dimensions with Gaussian interactions. We 
set the mean of the Gaussian interactions to Jo = —1 in 
order to suppress ferromagnetic ordering. 20 The effective- 
ness of this is shown in Fig. 1121 where the ferromagnetic 
susceptibility %fm is seen to be very small and does not 
exhibit measurable finite-size effects. However, if we take 
Jo = 0, which is also shown, xfm becomes large and size 
dependent, indicating that ferromagnetic ordering prob- 
ably takes place. The simulation parameters used for the 
ten-state Potts glass for Jo = — 1 are shown in Table [V] 
The second set of parameters for L = 6 and 8 is for addi- 
tional runs to probe the very-low-temperature behavior. 

A plot of £l (not divided by L) against temperature 
for Jo = —1 is shown in Fig. The correlation length 
shows no size dependence and stays very small, less than 
one lattice spacing, even at extremely low temperatures. 
We conclude, in agreement with Brangian et aZ.p2& that 
there is clearly no spin-glass transition in this model. 

Since the nearest-neighbor ten-state Potts glass does 
not even come close to having a transition down to T = 0, 
its behavior is very different from the mean-field (infinite- 
range) case, which has two transitions. Although the 
mean-field version provides a plausible model for super- 
cooled liquids as they are cooled down towards the glass 
transition, the nearest-neighbor model does not. It might 
therefore be interesting to investigate a Potts glass that 
has somewhat longer-range interactions^ to see if its be- 
havior corresponds more closely to that of supercooled 



liquids. 



V. CONCLUSIONS 

We show a summary of our results for the various mod- 
els in Table. IVll We find strong evidence for a finite tran- 
sition temperature T c in the three-state Potts glass in 
three dimensions, with universal critical exponents. This 
differs from several earlier studies which claimed that the 
lower critical dimension is equal or close to 3. In four 
dimensions, the three-state Potts glass has a transition 

TABLE VI: Summary of results for the four different p-state 
Potts glasses studied. The values for T c and v are obtained 
from analyzing data for (,l/L, which we argue give the most 
reliable results. The values for r\ are obtained from data for 
Xsg with only T c fixed to its value obtained from £l. 



V 


Disorder 


d 


T c 


V 


V 


3 


Gaussian 


3 


0.273(5) 


1.30(7) 


-0.03(2) 


3 


Bimodal (±J) 


3 


0.377(5) 


1.18(5) 


0.02(2) 


3 


Gaussian 


4 


0.536(3) 


0.81(2) 


-0.08(2) 


10 


Gaussian 


3 




No transition 





which is quite different from that proposed by Scheucher 
and Reger— — namely the discontinuous transition found 
in mean-field theory for p > 4. The ten-state Potts glass 
in three dimensions, which in its mean-field incarnation is 
often used to describe structural glasses, is very far from 
having a phase transition into a glassy phase at any tem- 
perature since its correlation length remains very small 
even down to extremely low temperatures, in agreement 
with the work of Brangian et alm& 
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